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ABSTRACT 



6 

in, 

c/3 , To investigate the physics of mass accretion onto weakly-magnetized neutron 

stars, 95 archival RXTE data sets of an atoll source 4U 1608-522, acquired over 
1996-2004 in so-called upper-banana state, were analyzed. The object meantime 
^ [ exhibited 3-30 keV luminosity in the range of < 10^^ — 4 x 10'^^ erg s~^, assum- 

Q I ing a distance of 3.6 kpc. The 3-30 keV PCA spectra, produced one from each 

25 ! dataset, were represented successfully with a combination of a soft and a hard 

|.^ I component, of which the presence was revealed in a model-independent manner 

O I by studying spectral variations among the observations. The soft component 

is expressed by so-called multi-color disk model with a temperature of ~ 1.8 
keV, and is attributed to the emission from an optically-thick standard accretion 
'k>( I disk. The hard component is a blackbody emission with a temperature of ~ 2.7 

H I keV, thought to be emitted from the neutron-star surface. As the total luminos- 

ity increases, a continuous decrease was observed in the ratio of the blackbody 
luminosity to that of the disk component. This property suggests that the mat- 
ter flowing through the accretion disk gradually becomes difficult to reach the 
neutron-star surface, presumably forming outflows driven by the increased radi- 
ation pressure. On time scales of hours to days, the overall source variability was 
found to be controlled by two independent variables; the mass accretion rate, 
and the innermost disk radius which changes both physically and artificially. 

Subject headings: stars: binaries: general — stars: individual(4U 1608-522) — 
stars: neutron — X-ray: stars 
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1. INTRODUCTION 

It has long been known that low-mass X-ray binaries (LMXBs), namely close binaries 
involving neutron stars (NSs) without strong magnetic fields, emit thermal type spectra 
when they are luminous (> 10^^ erg s~^). In the 1970's, their spectra were often modeled 
empirically in terms of simple thermal Bremsstrahlung, without much physical basis. As 
observations made progress, it gradually became clear that this simple empirical modeling is 
in fact inadequate, and at least two components are needed to describe wide-band (typically 
2-30 keV) LMXB spectra. Today, a canonical model (sometimes called "Eastern model") 
interprets the X-ray spectra of a luminous LMXB as originating from two major emission 
re gions; the surface of t he NS (or often called "boundary layer" ) and an accretion disk around 
it JMitsuda et al.lll984l ). 



The Eastern model was constructed based on Tenma observations of luminous LMXBs 



( Mitsuda et al.lll984lll989l : lMakishima et al.lll989l ). particularly considering intensity-correlated 
spectral changes. When an LMXB becomes brighter on a time scale of ~ 1000 s, its spec- 
trum hardens in energies from ~ 2 to ~ 10 keV, but the spectral shape in the hardest range 
(above 10 keV) stays constant. A difference between a pair of spectra with different inten- 
sities generally takes a form of a blackbody (BB) spectrum, of which the intensity changes 
but the shape (i.e., the temperature) is approximately constant. This BB component has a 
temperature of ~ 2 keV, which is close to the Eddington temperature, 2.0 keV, for a NS of 
~ 1.4 Mq mass (with Mq being the solar mass) and 10 km radius. In addition, the area of 
its emission region is inferred to be a fraction of the surface area of such a NS. Therefore, 
the BB component has been ascribed successfully to the NS surface emission. Indeed, every 
spectrum observed from these sources was reproduced by a sum of this BB component, and 
an additional soft component. The soft component was approximated first by a softer BB 
of ~ 1 keV temperature, but later, much better by a particular super position of BB spec- 



tra, called "disk blackbody" or "multi-color disk" (MCD) emission (JMitsuda et al.l Il984 



Makishima et al.l 119891) . as predicted by the optically-thick standard accretion disk model 
(jShakura fc Sunyaevlll973l ). 



While the Eastern model thus provides a promising ground to understand the physics 
of mass accretion in luminous LMXBs, their complex spectral and intensity variations have 
often been described in a purely empirical manner using "color-color" diagrams and other 
similar methods. As a result, such empirical classifications as "Z sources" and "atoll sources" 
have been created, together with various "branches" on these empirical diagrams (e.g., 
Hasinger &: van der Klislll989l ). These primitive descriptions are waiting to be replaced by 
more physically meaningful ones, using the Eastern model. 



After the launch in 1995, the Rossi X-Ray Timing Explorer {RXTE; iBradt et al.lll993l ) 
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has been observing many LMXBs for a huge number of times. With the largest effective area 
and the highest timing resolution ever achieved, RXTE is indeed very suitab le to the study 



of spe ctral yariations of LMXBs. On tira e scales of milliseconds to seconds, iGilfanov et al. 



( 120031 ) and iRevnivtsev fc Gilfanovl ( l2006l ) carried out such studies using the RXTE data, 
and revealed that the variation is carried by a rather hard BB-like emission component with 
a temperature of ~ 2-3 keV. These results reinforce the validity of the Eastern model, and 
encourage us to attempt its extensive application to the RXTE data of LMXBs. We thus 
hope to understand the spectral variations of LMXBs as a whole using a physical model 
(based on the Eastern description). 

As the first of a series of our planned publications, the present paper deals with so- 
called upper-banana state of 4U 1608-522, which is one of the most frequently observed 
atoll sources by RXTE. Section 2 describes 414 RXTE data sets of this LMXB, of which 95 
are selected for use in the present paper. In § 3.1, we select four representative spectra out 
of the 95 data sets, and analyze their intensity-correlated spectral changes to reconstruct the 
Eastern model. The derived model parameter values are examined in § 3.2, followed by § 3.3 
where the model is applied to all the 95 energy spectra and relations among the obtained 
physical parameters are studied. As a reconfirmation of our approach, we calculate in § 3.4 
the effective degrees of freedom involved in the variability, of which the results are physically 
interpreted in § 3.5. As discussed in § 4, the results give a full support to the Eastern model, 
and lead to a finding of mass outflows as the source gets luminous. 



OBSERVATIONS AND DATA REDUCTION 



Our target object 4U 1608-522 is an LMXB with a recurrent transient characteristic, 
exhibiting occasional outbursts. As the X-ray intensity evolves thr ough these outbursts, the 
sourc e is known to take all three spectral states of atoll sources (JHasinger fc van der Klis 
19891 ): namely islan d, lower-banana, and upper-banana states, in the increasing order of 
the source intensity (JMuno et al.l l2002l : iGierlihski fc Pond |2002aJ) . In the island and lower- 
banana states, a hard tail appears in the spectra presumably due to Com ptonization by hot 
electro ns, of which the temperature is of the order of several tens keV ( iGierlihski fc Done 
2002bl ). In the present paper, we analyze only the data in the upper-banana state, where the 
sp ectrum takes the ther mal-shape characteristic of luminous LMXBs. The original paper 
by iMitsuda et al.l (119841 ) also utilized 4U 1608-522, when the sourc e was presumably in th e 
upper-banana state. The source distance is assumed to be 3.6 kpc (JNakamura et al.lll989[ ). 



Although RXTE\ia,s observed many outbursts from 4U 1608-522, we limit the present 
study to those data sets which were obtained from the launch of RXTE (December 1995) 
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to the end of August 2004, when the brightest and the second brightest outbursts were 
recorded. During the period, 414 pointing observations of 4U 1608-522 were conducted 
with the Proportional Counter Array (PCA; iJahoda et aLlll996l ). Although some of the five 
Proportional Counter Units (PCUs) of the PCA were turned off in some observations, PCU2 
alone was always operational. In order to avoid systematic differences among different PCUs, 
we utilized the PCU2 data only, because the source in the upper-banana state is very bright 
(typically ~ 1000 counts s~^ per PCU) and hence statistical errors are usually negligible. 
We employed "Standard-2" data, of which the time resolution is 16 s, since they are most 
accurately calibrated and suitable for spectral analysis. The data were filtered in the standard 
manner for bright sources, and were corrected for dead times in the standard way. The PCA 
background was estimated for each dataset by PCABACKEST version 3.0. To remove Type 
I bursts, we excluded those time regions when the background-subtracted source count rates 
(per 16 s) are outside 50%-200% of the average over each observation. 

After the subtraction of the modeled background, we calculated color-color diagrams 
(CCDs) and hardness-intensity diagrams (HIDs) of 4U 1608-522 using data points averaged 
over 128 s. In this paper, the soft and hard colors refer to 6-10 keV vs. 2.5-6 keV source 
count ratios and those in 10-30 keV vs. 6-10 keV, respectively. The intensity used for 
HIDs refers to an energy band o f 2.5-30 keV in the unit of the Crab Nebula count rate. As 
described in lJahoda et al.l (|2006[ ). the overall PCA operation history is divided into 5 epochs 
according to differences in the high- voltage settings. As a result, CCDs and HIDs need to 
be produced separately for different epochs. Figure [T] representatively shows the obtained 
CCD and HIDs of 72 observations in the entire epoch 3. Here, we excluded faint data points 
of which the intensities are less than 10 counts s~^ PCU2~^ (^ 3 x 10^^ erg s~^ for 3.6 kpc). 

On the CCD of Figure [H the count rate increases from top through lower left to lower 
right, with each part corresponding to the island, lower-banana, and upper-banana states, 
respectively. On the HIDs, the data points form many vertical line-like features, or stripes. 
The source traces these stripes on a time scale of days, and moves horizontally on much longer 
time scales. The former causes the source motion within a single state on the CCD, while 
the latter leads to transitions among the three states. The hard color is higher in the island 
state, and decreases in the banana states. The present paper focuses on the upper-banana 
state, which is defined here as those datasets with the 3-30 keV source intensity exceeding 
0.4 Crab. We applied the same criteria to the data from the other four epochs. When 
summed over the 5 epochs, the selected datasets reach 95 in number, with their 3-30 keV 
luminosities spanning a range of (1- 4) xlO^'^ erg s~^. Some of t hese datasets ^yere alr eady 



analyzed with t he Eastern model by iGierlihski &: Dond (J2002bl ). iGilfanov et al.l (120031 ) and 
Lin et al.l f l2007l ). 
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For each dataset in the upper-banana state, we accumulated the Standard-2 data into 
a single 3-30 keV spectrum, with a typical exposure of several ks. Then, the modeled 
background mentioned above was subtracted. However, this method is known to causes 
slight (several percents) over-subtraction of the background |1|. This effect was compensated 
for by rescaling the background spectrum up to 10%, so that its count rate matches that 
of the on-source data in the hardest 60-100 keV energy range, where the signal count is 
considered negligible. Since the source count was detected significantly up to 30 keV, we 
hereafter analyze the total 95 energy spectra in the energy range of 3-30 keV. The PCU2 
response matrix was created for each observation using PCARSP version 10.1. In order to 
take into account calibration uncertainties, we add 1% and 7% systematic errors to each 
energy bin of the source and background spectra, respectively. 



. DATA ANALYSIS AND RESULTS 
3.1. Comparison of Energy Spectra 



In comparison with the previous work by iMitsuda et al.l ( 1l984l ) , we first examine rel- 
atively slow (> 10^ s) spectral variations along the upper-banana state of 4U 1608-522. 
For this purpose, we selected four representative spectra, hereafter Spec A through Spec 
D, obtained at the most luminous (several times 10^^ erg s~^) end (i.e., the right side) of 
the HIDs. Three of them (Spec A, Spec B and Spec C) come from the same vertical line, 
while the other is located on the different line. Their observation IDs, average intensities, 
and soft /hard colors are summarized in Table [H while their locations on the hard-HID are 
indicated in Figure [1] These spectra and their ratios to their average spectrum are shown 
in Figure [2l 

In Figure [21 the brightest (Spec A) and the second brightest (Spec B) spectra of the 
four are selected from the same vertical line in the HID. By comparing them, it is clear that 
the intensity in the hard energy band changes significantly, while that in the soft band is 
kept constant. Moreover, in energies above ^ 15 keV, the ir ratios saturates as a function 



of energy. As originally proposed by iMitsuda et al.1 (jl984l ). these results can be explained 



when the spectra consist of a variable hard component that dominates in > 15 keV energy 
band and a stable soft component, and the normalization (but not the shape) of the former 
changes between the two spectra. Indeed, the difference spectrum between them, shown in 
Figure [3] as Spec 1, is represented successfully by a BB model of which the temperature is 



^http://lheawww.gsfc.nasa.gov/users/craigin/pca-bkg-tdrift/ 
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kT^B ~ 2.5 keV. The obtained best-fit parameters are summarized in Table [21 Although the 



derived absorption column density, ^ 3x10 '^ cm , is significantly highe r than that reported 
previously (0.8 x 10^^ cm~^: lGrindlay Sz Lillerlll978l : lRutledge et al.lll999l ). the difference may 



be attributed to slight changes in the assumed stable soft component (as assessed later in this 
subsection), and /or those in /cTr r. These results reconfi rm the achievements with Tenma 
(JMitsuda et al.lll984l ) and Ginga ( jMakishima et al.l 119891 ). 



From the above results, we consider that the harder part of the spectrum is carried 
by a "hard component" which is app roximated by a ^Trr ^ 2.5 keV BB. Its variation (in 
normalization rather than in /cTrb; iMakishima et al.lll989[ ) causes significant changes in 
the hard color, but does not affect the intensity very much because the total source counts 
are dominated by softer photons. This explains the formation of the vertical stripes in 
the hard HID of Figure [H This hard component is most naturally attributed to optically- 
thick emission from the NS surface, because the measured temperature agrees with the local 
Eddington temperature of a 1.4 Mq NS, or ~ 2.0 keV. 

Contrary to the above case, the other two spectra, the faintest spectrum (Spec D) 
and the second faintest one (Spec C) which lie on a pair of stripes adjacent to each other, 
differ mainly in the soft energy band, with the hard energy end kept almost constant. This 
suggests that the source variation in this case is carried by a soft spectral component which is 
presumably identified as the stable component suggested by Spec A and Spec B. To constrain 
its spectral shape, we derived the difference spectrum (denoted Spec 2) between Spec C and 
Spec D, and show the results in Figure HI Clearly, this difference spectrum is much softer 
than the previous one (Spec 1), in agreement with the inference derived from the spectral 
ratios. This Spec 2 can be reproduced by either an MCD model {kTi^ ~ 1.7 keV) or a BB 
{kT-QB ~ 1.3 keV), but as shown in Table [21 the former is more successful and gives a more 
reasonable absorption. We hence consider that the "soft component" is approximated by an 



MCD model with kT„ 



1.7 keV. 



Since the above result supports the Eastern model, the four spectra must be reproduced 
by a linear combination (though with different weights) of the soft component expressed by a 
/cTin ~ 1.7 keV MCD, and the hard component identifiable with a /cTrr ~ 2.5 keV BB. Then, 
adding a narrow Gaussian component to represent an Fe-K line emission line at 6-7 keV, 
we fitted the four spectra with this composite model, denoted MCD+BB+Gau model. The 
results are shown in Figure [Hand Table [3l there, r^n is the innermost radius of the accretion 
disk for an assumed inclination angle 2 = 0°, and rRB is that of the BB model assuming 
a spherical emission region. Thus, all the four spectra have been fitted successfully by the 
model, and the obtained values of /cTin ~ 1.8 keV and /cTrr ~ 2.7 keV indeed agree with 
those obtained by the two difference spectra. These two temperatures are also consistent 
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with those found in pr evious works that employed the Eastern model (JMitsuda et al.lll984l : 
Makishima et al.lll989[ ). 



Can we explain these r esults with alternative modeling? For example, a model pro- 
posed by lWhite et al.l ( 1l988l ). often called "Western" model, assumed that Thomson opacity 
dominates free-free opacity in the accretion disk and the emission is significantly modified 
by electron scattering, while the emission from the NS surface are hardly modified and ob- 
served as a simple BB emission. Then, an unsaturated Comptonized emission is expected 
to dominate all over the energy band, while a soft BB model is added in the intermediate 
range. This model can actually explain the shape of Spec 2 in terms of variations of the 
soft BB component. However, Spec 1 is very difficult to explain with this model, because its 
shape is different from those of either component of the model; we would need extreme fine 
tunings among the parameters of the BB and Comptonized components. 

Returning to the Eastern model, a simple simulation was performed to examine the 
increased absorption issue in Figure [31 Employing the best fit model to Spec A, we created 
a fake spectrum to be called Spec A'. Yet another spectrum, named Spec B', was created, 
in which the BB normalization was reduced to 67% of that of Spec A', and the MCD 
normalization was set 5% higher than that in Spec A'. These are meant to emulate Spec B, 
within the allowed fit errors in Table Tabled Then, the difference between Spec A' and Spec 
B' was indeed fitted successfully by a single BB model with nearly the same BB temperature 
as assumed in the input, but with the absorption increased to ~ 3 x 10^^ cm~^. 



3.2. Absolute Values of the Physical Parameters 

The Eastern model allows us to assign clear physical meanings to absolute values of the 
physical parameters obtained by the MCD and BB components. For that purpose, however, 
we need to convert the measured color temperature kTx (X = "in" or "BB" ) to the effective 
temperature kT^, using so-called hardening factor k, as 

K,~^ kTx ■ 



kTf 



1 



Then, the true radius r|f is estimated from the measured value rx as 



„eff 



K rx . 



(2) 



The value of k is numerically es timated as 1.4-1.6 for Type I bursts (jEbisuzaki et al.lll984l ) 
and 1.7-2.0 for accretion disks (jShimura fc Takaharalll995l ). 



If we adopt k = 1.7, the observed soft- component parameters of kT^^ ~ 1.8 keV and 



4 km yields kT^ 



cff 



1.1 keV and r[ 



cflf 



12 km. Thus, the estimated true radius 



becomes larger than the NS radius of 10 km. Moreover, the value of r?f is close to the radius 
of last stable orbit in terms of general relativity, 3-Rs = 12.4 km, where Rs = 2GM/(? is the 
Schwarzschild radius, G is the gravitational constant, M is the NS mass, and c is the light 
speed. As to the BB component, /cTrr ^2.7 keV arid trr ^ 1.2 km yields kT^ ~ 1.8 keV 



and Tr^ ~ 2.7 km, assuming k ~ 1.5 (JLondon et al.lll986l : lEbisuzakil 119871 ) . The estimated 



kT^ is close to the local Eddington temperature (2.0 keV) of a 1.4 Mq NS, suggesting that 
the BB emission arises from a region on the NS surface. In addition, r^*^ is smaller than 
10 km, when assuming isotr o pic em ission from a spherical source. Therefore, as previously 
suggested by iMitsuda et all fll984[ ). the BB component can be regarded as being emitted 
from an equatorial zone of the NS, where the accretion disk contacts the surface. 



3.3. Relations among the Model Parameters 

From the above spectral analysis, it has been confirmed that the Eastern (MCD+BB+Gau) 
model successfully reproduces the selected four spectra in the upper-banana state, and yields 
physically reasonable interpretations. Given these, we applied the MCD+BB+Gau model 
to the 95 spectra in the upper-banana state, and obtained x^/d.o.f. of < 1.4 for all of them. 
Therefore, we regard the Eastern model as applicable to all the present data sets. 

Figure [6] summarizes the relation among the unabsorbed disk bolometric luminosity 
-^disk; the unabsorbed BB bolometric luminosity for isotropic emission Lrr, and their sum 
Ltot = -^disk + -^BB, as well as the temperatures and radii. Here, Ldisk is calculated assuming 
a face-on disk (i.e., the inclination angle i = 0°) as 



^disk 



/•oo 

2 X / 2TxraT{rYdr = Anrl^aT^ 



4 
in ) 



(3) 



where the first factor of 2 means emission from the two sides of the disk, a is the Stefan- 
Boltzmann constant, and T{r) is the disk temperature at the radius r. 

In Figures [6] (a), both Ldisk and Lrr are seen to increase as Ltot increases from 1 x 10^^ 
to 4 X 10^'^ erg s~^. However, a closer inspection reveals that Ldisk varies more steeply than 
Lrr. Actually, the data behavior in Figure E] (a) can be approximated as 



Lriisk oc L 



1.1 

tot' 



^BB 



OC L 



0.7 
tot' 



^BB 



OC L' 



0.6 

disk ■ 



(4) 



As a result, the ratio between Lrr and Ljisk de creases from 0.6 to . 4 as Lf^t increases by a 



facto r of 4, in agreement with previous reports (JGilfanov et al.ll2003l : iRevnivtsev &: Gilfanov 
20061 ). Similarly, Figure [6] (b) shows the luminosity dependences of the temperature and 
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radius parameters, derived in the same way as in § 3.1, and presented without applying any 
of the correction factors mentioned in § 3.2. The relations are approximated as 

kT,^ oc L°„f , ri, ex L°„f , /cTbb oc L°„f , tbb oc L^^^ . (5) 

Instead of Ltot-, we may utilize the mass accretion rate M, which can be estimated in 
the following way. Generally, the accreting matter releases energy at a rate of GMM/rin, as 
it flows through the disk down to the innermost radius Tin. According to the picture of the 
standard accretion disk, a half of this luminosity is radiated away in the accretion disk as 



-^disk = ^ , and the other half is stored in the Keplerian kinetic energy. Therefore, the 

■^ ^in 

mass accretion rate is estimated as 

^ = Gm''' ^ ^disknn oc rlj:^^. (6) 

This relation remains valid even if r^^ varies, as long as the disk stays in the standard state. 
Figure [7| shows the same physical parameters as presented in Figure El but this time, as a 
function of M estimated via equation ([6]), where the unit of M is arbitrary. 

In the case of a standard a ccretion disk, Lniay and fcTj^ are expected to be proportion to 



in ine case oi a siaiiaaiu a ccieiion uisk, i^hisV unu. iyi\n are expecieu lo ue pioporiion lo 
M an d M°'^^, respectively (e.g., iMakishima et al.lllQSGuEbisawa et al.lll993l : lTanaka fc Shibazaki 



19961 ). as long as nn kept constant. Indeed, in the range of M < 300 (corresponding to 



^tot 



< 



1.5 X 10^^ erg s ^), Figure [7] reveals tight scalings as 



Ldisk oc M0-92±0.07^ ^y.^ ^ ^0.21±0.05^ ^.^ ^ ^0.05±0.07 ^ (7) 

which agree, within er rors, with the pr edictions for standard disks. This result is consistent 



with that reported by iLin et al.l ( 120071 ). 



As the accretion rate M increases to > 300 in Figure [3, the disk parameters start 
deviating from the scalings of equation ([7]), and follow new relations approximated as 

Ldisk oc mO-80±o-03, kT,^ oc mO-ii±o-02^ ri„ oc mO-18±o.o3 _ ^8) 

Thus, the inner disk radius rjn apparently starts "retreating", with increasing fluctuations, 
but the increase in rin is compensated by the flattering in kT^^. As a result, /^disk starts to 
weakly saturate, instead of increasing in proportion to M. 

Turning to the BB parameters, we find /cTbb rather stable, with a weak positive depen- 
dence on M. On the contrary, tbb exhibits a large scatter, which causes similar scatter in 
Lbb- The M-dependence of these BB parameters can be approximated as 

Lbb oc M°•^ A;Tbb oc mO-07±o.oi^ ^^^ ^ ^o.is _ ^g) 
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This BB behavior is the same as previously reported (JMitsuda et al.lll984t iGilfanov et al. 
20031 : iRevnivtsev fc Gilfanovll2006[ ). We also find that /cTbb is always higher than kT[^, and 
r-Q-Q is smaller than rjn. These relations are consistent with the basic assumptions of the 
Eastern model; the BB emission comes from the NS surface, while the MOD component 
from the surrounding accretion disk. 



3.4. Effective Degrees of Freedom Causing Spectral Variability 



We have so far described the energy spectra and their variations in terms of the four 
quantities, fcTjn, r^n, /cTbb, and tbb- However, it is not yet clear how many of them can be 
regarded as independent variables (with the rest depending on them). In other words, we 
need to know how many degrees of freedom are involved in the observed spectral variations 
in the upper-banana st ate. This can be done by applying "fractal dimension analysis" (e.g., 
Matsumoto et al.ll2005l ) to the four model parameters. 



Yl — n.-Lin, 



For this purpose, let us define a 4-dimensional vector space spanned by the four variables, 
: fcTin, Y2 = Tin, Fs = /cTbb, and I4 = tbe, wherein each data set is represented by a 
vector 

Vlt) = {^(z), Y^it), Ysii), Y,it)} 

with i = 1,2, ..,95 denoting the data number. Let us also define the average vector {V) = 
E^iil^(z)/95 = {{Yi),{Y2),{Y3),{Yi)}, and the zero-mean vectors, v{i) = V^i) - (V) = 
{yjii) , y2{i) , y3{i) , y4{i)} with yk{i) = Yk{i) - (Yk). Then, a "distance" D{{) of the vector 
v{i), measured from the origin, is calculated as 



DH) 



1I/2 



^{yk{i)/(rky 



k=l 



l,2,---,95), 



(10) 



1/2 



where ak = [X^j=i 2/fc (0^/94] ' is the standard deviation in yk{i) around the origin [or in 
Yfe(z) around its mean (Yk)]. Finally, we calculate the number A^(< D) of those data points 
of which the distance D{i) is less than a given value D. 



Figure [8] shows the normalized data point number A^(< D)/95 from equation (fTOl) as a 
function of D, over the range of A^/95 = 0.2 — 0.8 (or A^ = 23 — 76). If the variations are 
controlled by n (1 < n < 4) independent parameters, we expect the vectors v{l),v{l), ...v{95) 
to form a n-dimensional subspace in the vector space, so that A^(< D) should increase as 
oc D^. Indeed, Figure [8] reveals a tight power-law relation as 



N oc Di-8±o-i 



'11^ 
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Since this result is consistent with ra = 2, we infer that the spectral behavior in the upper- 
banana state of 4U 1608-522 has effectively two degrees of freedom. 



3.5. Fluctuations Independent of the Total Luminosity 

Of the two independent variables describing the spectral variability (§ 13.41) . one is obvi- 
ously the total luminosity Ltot, or nearly equivalently, the mass accretion rate M . Actually 
in Figure 6 (or Figure 7), the four quantities are all observed to depend primarily on Ltot (or 
M). However, they show significant scatters around the Ltot-dependent correlation, so that 
none of them can be regarded as a single- valued function of Ltot (or M). This is consistent 
with the presence of the second degree of freedom revealed in § 13.41 In order to identify what 
is causing this extra freedom, we remove the Ltot-dependence from the behavior of the four 
model parameters, and study the residual variations. 

To eliminate the major Ltot-dependence from the physical parameters, we calculated 
their de-trended counterparts Z\i) {Z = Ldisk, Lbb, fc^in,nn, ^^bb and tbb), employing 
equations (jlj) and ([5]), as 

LU(^) oc Ldisk(^)/Ltot(0'-\ ^bb(0 ^ LBB(0/^tot(^)°-', (12) 

kTf^i^ oc A;^i„(^)/Ltot(^)°■'^ C{^) oc rU^/LU^f-'^ (13) 

fcr^B(OocA;TBB(0/Ltot(0°-'°, and r'^j,{i) (x rBB{^) / Ltot{^f■'' ■ (14) 

(As the de-trending parameter, we chose Ltot rather than M, because the latter is less directly 
estimated from the data.) Absolute values of the de-trended parameters are not meaningful, 
and we hereafter consider their relative values only. Since r[^ shows the largest variation 
among them, we plot in Figure [9]the de-trended parameters as a function of r{^. Table H] also 
summarizes those of Spec A to D in § 13. 1[ From the behavior of the de-trended luminosities, 
the data points can be divided into two branches; one has almost constant luminosities as r[^ 
varies, while the other is characterized by significant r(^-dependent variations both in L^j^j^ 
and Lgg. Hereafter, the two branches are denoted as "constant-luminosity branch (CLE)" 
and "variable-luminosity branch (VLB)", respectively. The two branches may connect at 
r{^ ~ 0.85, rather than behaving independently from each other. 

In the CLE, the de-trended EE luminosity stays nearly constant at Lgg ~ 0.8; so are 
the two EE parameters, kT^^ ~ 1.6 and r^g ~ 0.9. Similarly, the de-trended disk luminosity 
remains at L^^gj^ ~ 0.9. As a result, the CLE data points are distributed in Figure E] (a) 
along the major correlation trends. Nevertheless, r'^^ varies by ~ ±20%, accompanied by 
a clear decrease in the de-trended disk temperature as kT(^^ oc rin~°'^. (This scaling is a 
natural consequence of the constant L^j^j^ and the relation of L'^^^^^ oc r{^^^T(^'^.) There are two 
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possibilities to explain this r[^ behavior. One is real changes in rjn, and the other is those 
in the color hardening factor k of the disk. Since L'^^^^^ is kept constant in the CLB, the 
latter case may be more likely. Specifically, the observed behavior will be explained if some 
unspecified mechanisms (possibly related to vertical structure changes in the disk) caused 
K to increased by ~ 10% while somehow keeping the effective temperature of the disk and 
its true radius constant; the observed color temperature kT(^ would then increase by ~ 10% 
according to equation ([T]), and the apparent disk radius r(^ would decrease by ~ 20% after 
eqation (E]), thus reproducing the observaton. 

In the VLB, r[^ varies over a relatively small range (~ ±10%), whereas the two lumi- 
nosities both vary significantly in an anti-correlated way; L'^i^k decreases from 1.9 to 1.5 as 
oc r[j^~^, while Lgg increases from 0.8 to 1.2 as oc r[^^. This complementary behavior between 
the MCD and BE components appears in Figure M as the large fluctuations away from the 
major trend, where the data points reach a level of L^isk ~ -^bb ~ O.SLtot- (This behavior 
cannot be a fitting artifact caused by strong couplings between the two spectral component, 
since the nominal fitting errors are not particularly larger in the VLB, and are generally 
smaller than/comparable to the size of the plotting symbols in Figure [H] and Figure [71 ) 
Since /cTgg is almost constant (or only slightly decreasing) like in the CLB, the increase 
of Lgg is mainly caused by that of r^g from 0.9 to 1.2 as oc r{^^. We again observe kT(^ 
to decrease like in the CLB, but with a different r[j^- dependence whic h is approximated as 



oc r- . This is exactly the behavior of a standard accretion disk (iShakura fc Sunyaev 



19731 ) when its radius varies under a constant M, because we then expect M oc rf^T^ (equa- 
tion ([H])). This also explains the observed disk luminosity behavior, as -^^isk °*^ '^fn^in^ oc r[^^^. 
Obviously, the increase in the disk radius under a constant M predicts that a larger fraction 
of the overall gravitational energy release should be emitted by the BB component; in fact, 
along the full VLB in Figure Et^b), L'^isk is seen to decrease by ~ 0.3 in the employed arbitrary 
unit, while Lgg to increase by ~ 0.3, thus conserving the total luminosity. In other words, 
the VLB is characterized by actual changes in the innermost disk radius. 

Based on these correlation analyses among the de-trended spectral parameters, we con- 
clude that the second independent variable (besides M; § 13. 4p causing the variability in the 
upper-banana state can be ascribed to sporadic changes in r[^ (or in ri^). However, this r[^ 
variability is likely to be further subdivided into two different mechanisms under a constant 
M; one is an apparent effect due to variations in k of the disk, and the other is real changes. 

As shown in Figure [9] and Table HJ the L'^^^i^ parameters of Spec B, C and D locate on 
the CLB line, while that of Spec A does on the VLB one. Therefore, we consider that the 
hard/soft-band differences of Spec A/B and C/D observed in Figure |2] are mainly attributed 
to the VLB and CLB variation, respectively. On the other hand, the de-trended parameters 
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of Spec B and C are similar, and the variability between them observed over the entire 3-30 
keV band is thought to be dominated by the M change. 



4. DISCUSSION AND CONCLUSION 

4.1. The Eastern Model 

Through the analysis of the i?XT£' spectra of the atoll source 4U 1608-522 in its upper- 
banana state, we have confirmed that both time-averaged and difference energy spectra 
can be reproduced successfully by the Eastern (MCD-|-BB+Gau) model. Considering the 
hardening factor, the effective temperature and radius of the MCD component were obtained 
as kT^^ ~ 1.1 keV and r?f ~ 12 km, respectively. The radius is larger than the representative 
NS radius, 10 km, and is consistent with the last stable orbit, 3-Rs = 12.4 km, allowed by 
general relativity. This result agrees with the picture of a standard accretion disk which is 
formed around a NS. The BB parameters were found as kT^ ~ 1.8 keV and r*^^ ~ 2.7 km, 
assuming isotropic emission from a spherical source. The temperature is close to the local 
Eddington temperature (2.0 keV) at the NS surface, and the radius is smaller than 10 km. 
Then, the BB component can be regarded as being emitted from an equatorial zone of the 



NS, as previously suggested by iMitsuda et all (Il984[ ) 



4.2. Saturation of the Blackbody Luminosity 

According to the picture of standard accretion disks, a half of the released gravitational 
energy is radiated from the accretion disk (//disk)? and the other half is stored in the Keplerian 
kinetic energy and then emitted (Lbb) when the matter settles onto the NS surface. Then, 
we expect Lbb to be proportion to Ldisk, and hence the LBB/-^disk ratio to be constant. 
Indeed, Ldisk was found to increase almost linearly with the total luminosity Ltot (figure 6a). 
However, the same figure reveals that Lbb increases less steeply, making the LBB/Ldisk ratio 
decrease from 0.6 to 0.4 as L^t increases from 1 x 10'^^ to 4 x 10'^^ erg s~^. 

We may think of two possibilities to explain the observed relative decrease of Lbb- One 
is that the accreting matter reaches the NS surface and emits Lbb which is equivalent to the 
Keplerian energy (i.e. oc Ldisk), but we cannot observe all the emission because its wavelength 
shifts outside the PCA energy band (3-30 keV), or the geometrical angle of the emission 
moves from our line of sight. However, we have not observed any hint of extra emission in 
the softest or hardest energy ends of the PCA spectra. In addition, other sources, which are 
thought to have different inclination angles, are also reported to show the same behavior of 
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the decreasing LBB/-Z^disk ratio fJRevnivtsev fc Gilfanovll2006l ). Therefore, this possibihty is 
unhkely. 



The other possibihty is that the accreting matter flows aU the way through the disk 
down to its inner radius that is close to the NS surface, but its progressively larger fraction 
fails to accrete onto the NS surface. If the accretion-failed matter stayed around the NS away 
from the surface (e.g., an expansion of the boundary layer), it would accumulate to become 
optically thick, eventually producing detectable emission. This would lead to a decrease in 
/cTbb, because the emission radius should increase. As shown in Figure Eland equation ([S]), 
however, this is not the case; the kT^B value stays almost constant (or rather increases) as the 
total luminosity increases. Then, a fraction of the matter must be outflowing and escaping 
from the system without releasing its kinetic energy as emission. Although the observed 
Ltot ~ 4 X 10^'' erg s"^ is only ~ 20% of the Eddington luminosity for a 1.4 Mq NS, the 
observed value of kT^ ~ 1.8 keV is already close to the spherical Eddington temperature at 
10 km (2.0 keV). Considering further the non-spherical geometry of the disk, and the disk 
radius which is possibly larger than 10 km, the putative outflow is likely to be driven by 
increased radiation pressure in the innermost disk region. The presence of such outflows is 
indeed suggested by detections of broad absorption features from LMXBs (Schulz & Brandt 
2002; Ueda et al. 2004), black hole binaries (Kotani et al. 2000; Yamaoka et al. 2001; Kubota 
et al. 2007), and active galactic nuclei (Pounds et al. 2003ab). 

Assuming that the decrease in LBB/-^disk from 0.6 to 0.4 is due to outflows, about 
(0.6 — 0.4)/0.6 ~ 30% of the accreting matter is estimated to escape from the system at 
a typical luminosity of ~ 4 x 10^^ erg s~^. Toward lower luminosities, the ratio appears 
to converge to ~ 0.6, rather than 1.0 which would be theoretically expected. This may be 
attrib uted to effects due, e . g., to the system inclination, general relativity, and NS rotation 



Sunyaev fc Shakural (1l986l )). Detailed theoretical discussion on this point is beyond 



(e.g., 

the scope of the present paper. 

As indicated by equation ([8]), the innermost disk radius r[^ is observed to increase 
gradually when the mass accretion rate increases to> 300 in Figure [3, or equivalently, when 
the total luminosity becomes > 1.5 x 10'^^ erg s^^ (corresponding to > 7% of the Eddington 
luminosity). This is unlikely to be an apparent effect caused, e.g., by chang es in the hardening 



facto r , since this quant ity would increase as the luminosity increases ( iGierlihski fc Done 



20031 : iDavis et al.ll2005l ). and would hence make r^^ smaller. Therefore, the increase in rjj 



toward higher accretion rates is considered to represent a real retreat of the innermos t disk 



radius. This interpretation may agree with the results by iPopham fc SunyaevI (120011 ) . who 
numerically calculated the physical condition of the boundary layer and showed that the disk 
inner edge retreats back by the increased radiation pressure when the luminosity approaches 
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the Eddington limit. 



4.3. Secondary Spectral Variations 

The result of the fractal dimension analysis shows that the source variations are con- 
trolled by two independent variables. One of them is obviously the mass accretion rate (or 
nearly equivalently the total luminosity), which produces the major trend in Figure [7] as dis- 
cussed so far. The other variable is thought to be random changes in the inner disk radius, as 
represented by the behavior of r[^ in the de-trended analysis performed in § 13.51 As discussed 
therein, the variations of r[^ may be caused by two different mechanisms, corresponding to 
the CLB and VLB. These two branches are considered to degenerate in the fractal dimension 
analysis, and difficult to identify as two independent degrees of freedom. 

In the CLB, the change in r[^ seems to be apparent rather than real, and is considered to 
reflect fluctuations in the hardening factor of the disk As a result, the BB parameters remain 
unchanged, while the MCD parameters vary as kTl^ oc r[^~°'^, keeping the disk luminosity 
constant (L^j^j^ oc r[^). Such variations in the hardening factor may in turn stem, e.g., from 
fluctuations of the disk height, because this would affect the temperature gradient and/or 
radiative transfer in the vertical direction of the disk. 

In the VLB, the change in r[^ is considered real (§ l3.5p . because we observe a scaling of 
kTl^ oc r[^~^-'^^ which agrees with the behavior of a standard accretion disk under a constant 
M and a constant hardening factor. In this branch, the disk is considered to randomly 
retreat back from its "home position", which is located in Figure IH] around r[^ ~ 0.85 where 
the CLB and VLB meet. When this takes place, the accreting matter located between the 
retreated and the home-position radii will fall onto the NS surface without radiating the 
disk emission, and the gravitational energy it acquires will be released solely as the radiation 
from the NS surface. Then, L^isk will decrease while Lbb will increase. This can naturally 
explain the anti correlation between the MCD and BB luminosities observed in the VLB 
(Figure ED. 



4.4. Conclusion 

The results of our study can be summarized in the following four points. 

1. The physical picture provided by the Eastern model can consistently explain the be- 
havior of 4U 1608-522 in the upper-banana state. 
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2. When the mass accretion rate increases, some fraction of the matter flowing through 
the accretion disk starts outflowing, presumably due to an increased radiation pressure. 

3. Even at the same accretion rate, the spectral parameter changes randomly on time 
scales of hours to days, mainly reflecting changes in the innermost disk radius. 

4. The change in the disk inner radius is a mixture of real and apparent effects, the 
former due to sporadic retreat of the disk while the latter due to fluctuations in the 
color hardening factor. 

Facilities: RXTE (PCA). 
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Table 1: Four spectra analyzed circumstantially in the upper-banana state. 
Spectrum Observation Count Rate Color 

ID [Crab] Soft'' Hard^ 



Spec A 30062-03-01-02 1.07 0.43 0.29 

SpecB 30062-03-01-01 TOS 041 0.26 

Spec C 30062-03-01-00 (195 040 0.25 

SpecD 30062-03-01-04 (TtS 042 0.27 

'^ The soft color refers to the source count ratio of 6-10 keV / 2.5-6 keV. 
^ The same as the soft color but the ratio of 10-30 keV / 6-10 keV counts. 



This preprint was prepared with the AAS IATJt;X macros v5.2. 
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Table 2: Model fit results to the difference spectra calculated from those in Tabled]^ 
Spectrum Model A^h ^ kT^B ' kT,^ " xVd.oi~ 

Spec 1 BB 27+n 2.5 ±0.1 ~ 20/58 

(A-B) 

Spec 2 BB <0.3 1.27 ±0.02 ■■■ 47/58 

(C-D) MCD 1.8 ±1.2 •■■ 1.7 ±0.1 29/58 

^ All the errors are single-parameter 90% confidence limits. 
^ The absorption column density in the unit of 10^^ cm~^. 
'^ Temperatures are all in the unit of keV. 



Table 3: Fitting results of the four original spectra in Table [T] with the MCD±BB±Gau 
model. ^ 



Spectrum 


MCD 




BB 




Gau^ 






fciin 


'in 


fcTBB^ 


r-BB "" 


Cen. Ene.^ Norm§ 


xVd.o.f. 


Spec A 


1 O/I+0.03 


4 1+0.1 

^•J^-0.2 


2.67 ±0.04 


1.5 ±0.1 


6.6 ±0.2 1.9l[]j 


29/54 


Spec B 


1.88 ±0.04 


A n+0.1 


9 78+0-1 
^•'^-0.08 


1.1 ±0.1 


6.6 ±0.2 2.0 ±0.4 


41/54 


SpecC 


1.79 ±0.04 


A O+0.2 


2.64 ±0.08 


1.2 ±0.1 


6.6l|];^ 1.7 ±0.4 


32/54 


SpecD 


1.84 ±0.06 


q c:+0.3 

0.d_Q2 


2.7 ±0.1 


1 1+°-2 


6.6 ±0.2 1.11°;^ 


42/54 



'"^ All the errors are single-parameter 90% confidence limits. The absorption column density is fixed at 

0.8 X 10^2 cm-2. 

"" The width of the Gaussian component is not well determined and constraint to be < 0.2 kcV. 

'^ Temperatures are all in the unit of keV. 

'^ The innermost radius in the unit of \Jcosi km. where i is the inclination angle. 

" The radius of the BB emission in the unit of km, assuming the emission comes from the isotropic 

sphere region. 

^ The center energy of the Gaussian model, which is limited in the range of 6.4-6.9 kcV. 

s The normalization of the Gaussian model in the unit of 10^^ photons cm^^ s^^. 
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Table 4: de-trended parameters of Spec A to D after eliminating the Ltot dependence in 
Figure E"'^ 





-^disk 


^BB 


kTL 


r'. 

in 


^^BB 


'^BB 


Spec A 


1.69 


1.11 


1.43 


0.97 


1.58 


1.15 


Spec B 


1.96 


0.75 


1.49 


0.96 


1.65 


0.86 


Spec C 


1.93 


0.80 


1.44 


1.04 


1.58 


0.97 


Spec D 


1.90 


0.85 


1.53 


0.88 


1.64 


0.92 



^ See equation (jlip for the Lt^t-dependence eliminated here. 

^ Absolute values of the de-trended parameters are not meaningful. 
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Fig. 1.— The CCD (panel a), soft-HID (panel b), and hard-HID (panel c) of 4U 1608-522 
in epoch 3, with the time bin of 128 s. The utilized energy bands are described in the text. 
Four arrows indicate the datasets analyzed in 5 3.1 (Tabled]). 



-22- 







^ ■ ' ' ■ — 1 


A : 


I 


r 


^^ 


B ■ 
C 1 




' 


^%. 


D : 






■^^ 


■ 


"^ o 




^• 




> ^ 


■ 


■5 




^ 






^- : 


Counts s" 
01 0.1 1 


r 

(a) 


" 





1.7 ■ 



1.3 ■ 



0.9 
0.8 
0.7 
0.6 • 



0.5 



10 
Energy (keV) 



20 




^. _ _ '" "?'*'"i*;;^i* • 



D 



(b) 



TO 

Energy (keV) 



20 



Fig. 2. — Four energy spectra (panel a) shown in Table [T] and Figure [T] (c), and their 
intensity ratios to their average (panel b). The energy spectra are shown without removing 
the detector response. 
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Fig. 3. — The difference spectrum (Spec 1) between the brightest (Spec A) and the second 
brightest (Spec B) spectra in Figure [2l It is successfully represented by an absorbed BB 
model with a temperature /cTbb ~ 2.5 keV (lower histogram). The fit residuals are shown 
in the bottom panels. If the absorption is fixed at 0.8 x 10^^ cm~^, the model over-predicts 
the data below 7 keV (upper histogram). 
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Fig. 4. — The same as Figure |3l but for the difference (Spec 2) between the faintest spectrum 
(Spec D) and the second faintest one (Spec C) in Figure [21 It is reproduced by either a BB 
model (panel a, /cTbb ~ 1.3 keV) or an MCD model (panel b, /cTjn ~ 1.7 keV), with the 
absorption column density left free. 
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Fig. 5. — The same four spectra as presented in Figure [2] (a), fitted by an MCD+BB+Gau 
model with the absorption column density fixed at 0.8 x 10^^ cm~^. 
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Fig. 6. — (a) The relations of the MCD luminosity (Ldisk) and the BB luminosity (Lbb) as 
a function of L^ot = -^disk + -^bb in the upper-banana state. Each data point represents one 
of the 95 observations. The three dashed lines show O.TLtot? 0.5Ltot, and O.SLtot, from top 
to bottom, (b) The same as (a) but for the temperatures and radii of the MCD and BB 
models. 
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Fig. 7. — The same as Figure El but sliown as a function of the estimated mass accretion rate 
M. The three dashed hues illustrate dependences as Ldisk oc M, Hn oc M°, and kT\^ oc M°-^^. 
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Fig. 8. — Results of the fractal dimension analysis of the four model parameters {kT^a.-, ^h 
kTBB, and rBe) over the 95 data sets. See text for details. 
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Fig. 9. — The same as Figure Ej but for the physical parameters from which the Ltot 
dependence was removed via equations (12)-(fT4|). They are plotted as a function of the 
de-trended innermost radius of the accretion disk r[^. Dashed lines indicate the relations of 
r[^^, with p = —1.0, —0.75, —0.5, 1.0 and 2.0. The L'^is]^ parameters of Spec A to D in TableH] 
are also indicated. 



